Non-equilibrium dynamics of Andreev states in the Kondo regime 
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The transport properties of a quantum dot coupled to superconducting leads are analyzed. It 
is shown that the quasiparticle current in the Kondo regime is determined by the non-equilibrium 
dynamics of subgap states (Andreev states) under an applied voltage. The current at low bias is 
suppressed exponentially for decreasing Kondo temperature in agreement with recent experiments. 
We also predict novel interference effects due to multiple Landau-Zener transitions between Andreev 
states. 
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The interplay between Kondo effect and superconduc- 
tivity can lead to a rich phenomenology like the exotic 
behavior of heavy- fermion compounds [l| . Artificial de- 
vices like quantum dots (QD's) now open the possibility 
of exploring this interplay in well characterized trans- 
port experiments Very recently, Buitelaar et al. Q 
have achieved a physical realization of a S-QD-S system 
using carbon nanotubes coupled to Al/Au leads. Their 
results, showing an enhancement of the conductance in 
the Kondo regime with respect to the normal case, have 
not hitherto received a theoretical explanation. 

Most theoretical efforts on S-QD-S systems have been 
so far devoted to analyze the Josephson transport in equi- 
librium conditions ^ .5]. This system exhibits a transi- 
tion from a Kondo-like to a localized moment regime as 
a function of Tk /A, where Tk is the Kondo temperature 
in the normal state and A is the superconducting energy 
gap. Associated with this transition the current-phase 
relation undergoes a reversal of sign usually known as 
TT-junction behavior ^|^. 

A remarkable feature of this system, which is common 
to all coherent nanostructures coupled to superconduct- 
ing electrodes, is the appearance of subgap states (usu- 
ally known as Andreev states (AS's)) which give the main 
contribution to the Josephson current . These current 
carrying bound states arise from an infinite series of An- 
dreev refiections at the superconducting electrodes. In 
Ref. it was shown that electron correlations in a S- 
QD-S in the Kondo regime give rise to a renormalization 
of the AS's which is consistent with the renormalization 
of the Fermi liquid parameters in the normal state. 

In order to obtain theoretical predictions that could be 
compared with actual transport experiments one should 
address the situation in which a finite bias voltage is 
applied. This is a challenging theoretical problem in 
which the current is due to multiple Andreev reflection 
(MAR) processes involving strongly correlated electrons. 
Although some numerical calculations have already been 
presented 0, a detailed understanding of the transport 
properties in the voltage biased situation is still lacking. 
The problem is particularly involved in the low bias limit 
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FIG. 1: Schematic representation of the Andreev states dy- 
namics in two different situations: a) quantum point contact 
and b) S-QD-S in the Kondo regime. The dashed lines cor- 
respond to the states for perfect normal transmission. The 
lower and upper striped regions represent the states in the 
continuous spectrum. The straight arrows indicate Landau- 
Zener transitions between AS's, while the curved arrows rep- 
resent inelastic transitions between the AS's and the contin- 
uous spectrum (see text). 



in which the current is due to Andreev processes whose 
order n ~ int(2A/ev) becomes arbitrarily large as the 
voltage V tends to zero. The aim of this letter is to show 
that these properties can be understood in terms of the 
dynamics of the Andreev states under an applied bias 
voltage. 

In order to give a qualitative explanation of this dy- 
namics let us first analyze the behavior of the AS's in dif- 
ferent situations, as illustrated in Fig. 1. In the absence 
of a QD, a single channel S-S contact with normal trans- 
mission T exhibits AS's located at ±Ay^ 1 — r sin^ (j)/2 
0, where is the phase drop across the contact. In 
equilibrium and at low temperatures only the lower state 
is occupied. Under a small bias voltage this state evolves 
adiabatically except for Landau-Zener (LZ) transitions 
to the upper state with probability p — exp(— 7rAr/ew), 
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where r = 1 — r This mechanism produces a net 
transfer of quasiparticles from the occupied to the empty 
states above the superconducting gap accounting for the 
low bias dc current in the system, as represented by the 
arrows in Fig. 1 a. As shown below, in a S-QD-S in the 
Kondo regime the situation is qualitatively different due 
to the fact that the AS's detach from the continuous part 
of the spectrum (see Fig. lb). As a consequence the con- 
version of bound state electrons into free quasiparticles 
now requires an inelastic transition across the gap be- 
tween the AS and the edge of the continuous spectrum. 
This feature is the main ingredient which determines the 
behavior of the current at low voltages. In particular 
we will show that there is an exponential suppression of 
the current as a function of the parameter Tft-/A. The 
detachment of the AS's from the continuous spectrum 
also gives rise to novel interference effects due to multi- 
ple LZ transitions between the bound states which could 
be observed directly as an oscillatory behavior of the dc 
current. 

To describe a S-QD-S we start from the usual sin- 
gle level Anderson model, thus assuming a sufficiently 
large level separation. The model Hamiltonian is H — 
Hl + Hr + Ht + Hd, where H^ ji describe the 
left and right leads as BCS superconductors, Ht = 
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- h.c, is the term coupling the dot 

to the leads, d], and ^ being creation operators for 
electrons in the dot and in the leads respectively. Hd = 
Eo- ^of>'a-'rUii'\iii is the uncoupled dot Hamiltonian char- 
acterized by the dot level position eo and the Coulomb 
interaction U . It is assumed that the coupling of the 
dot to the leads in the normal state can be described by 
energy independent tunneling rates F^ = F^^ = F. As 
discussed in Ref. ^3 the bias voltage can be introduced 
as a time-dependent phase factor in the hopping term. 

A crucial point in the analysis of this model is how to 
deal with the electron correlation effects. In Ref. a 
perturbative approach was used to study the zero- voltage 
case. It was found that the AS's in the Kondo regime 
{Tk > A) have essentially the same phase dependence 
as in the non-interacting case but with a reduced ampli- 
tude. This is consistent with a Fermi liquid description 
of the normal state with renormalized parameters eg and 
F*, where F* ~ fixes the width of the Kondo res- 

onance. Although different approximations differ in the 
way to relate these quantities to the bare parameters, the 
description in terms of the renormalized ones can be con- 
sidered as universal ^he opposite hmit {Tk < A), 
the system is well described within the Hartree-Fock ap- 
proximation 0,0- 

In the present work we shall be concerned with the 
Kondo regime. At zero bias the AS's are determined by 
the poles of the dot retarded Green function, which in 
the Nambu space and close to the Fermi energy adopts 
the form 
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— {to + iO+)^ correspond to the dimensionless 
Green functions of the uncoupled electrodes. In the 
electron- hole symmetric case {eq = 0) the AS's are then 
determined by the equation 
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The solutions of Eq. (0) for F* > A and |0| < 1 are 
well approximated by cOs = ±A cos 0/2 with A = A[l — 
2(A/2F*)^]. Thus, as commented above, the AS's for 
the S-QD-S are detached from the continuous spectrum 
for finite values of T^/A, with a gap given by A — A 
(see Fig. 1 b). For decreasing T^/A the AS's begin to 
deviate from this simple phase-dependence 

In a non-symmetric situation {eq ^ 0) the two ballis- 
tic states are coupled and there appear an upper and a 
lower bound state with an internal gap of the order of 
2AV1 - T, where t = 1/ [l-l- (e;5/2F*)2] is the normal 
transmission at the Fermi energy. This is similar to the 
behavior of the AS's in a point contact with finite re- 
flection probability, except for the detachment from the 
continuous spectrum (see Fig. 1 b). 

In order to analyze the voltage biased situation we shall 
assume that the description in terms of the renormalized 
parameters Eq and F* is still valid in the voltage range 
{ev <C A) in which we are interested. Notice that in the 
Kondo regime this corresponds to ev <C Tr- and thus one 
can safely assume that the coherence associated with the 
Kondo effect will not be destroyed [l^ . 

For calculating the current we use a combination of 
Keldysh and Nambu formalisms following Refs. pTij |. 
Explicit extensions of this approach to the case of in- 
teracting quantum dots are discussed in Refs. 0, Il3| . 
This procedure allows to evaluate the current in terms of 
Green functions as I{t, v) = ^n{v) exp [m0(t)], where 
(j>{t) = 2evt/h -1-00- In the present case we shall be in- 
terested in the dc component Iq {v) . 

We shall analyze first the electron-hole symmetric case 
(e^ = 0). The behavior of the dc current at low bias is 
shown in Fig. 2. When F* ^ A the IV characteristics of 
a ballistic quantum point contact is gradually recovered, 
exhibiting a divergent low bias conductance with a satu- 
ration of the current &i v = towards the value 4eA//i 
mi^. As F*/A is reduced the low bias conductance is 
gradually suppressed. When F* A oscillations with a 
period A/n, where n is an integer, starts to be observable 
in the IV characteristic. These oscillations arise from the 
opening of new Andreev channels. As discussed in Ref. 
MAR processes of odd (even) order are enhanced 
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FIG. 2; Low bias IV characteristics of a S-QD-S in the Kendo 
regime for perfect normal transmission and different values of 
T'/A = 1, 2, 3, 4, 5 and 10 from bottom to top. The full 
lines correspond to the numerical results, while the dashed 
lines correspond to the approximation given by Eq. 



(suppressed) due to the presence of the dot level at the 
center of the gap. 

These results can be understood in more detail by an- 
alyzing the dynamics of the AS's under a finite bias volt- 
age. In order to establish a dc current the quasiparti- 
cles should be excited from the AS's to the continuum 
states. A simple diagonalization of Eq. for Eq = 
allows to identify = uj\/ A'^ — uj^/2r* as a complex 
self-energy for the ballistic states at ±A cos 0/2. The 
real part of within the gap is responsible for the 

renormalization of the bound states for finite F* (see Eq. 
(0)). When a small bias voltage is applied the AS's tend 
to follow adiabatically the evolution of the phase. The 
transition rate into the continuum is determined by the 
imaginary part of the self-energy and can be calculated 
using standard time-dependent perturbation theory. The 
probability for this transition is then given by 
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Using the approximate form of the AS, — A cos 0/2 
(valid for large r*/A) and in the stationary limit {t — > oo) 
one obtains the rate for this transition, jt ~ X]n>oT"' 
with 



7„ = ^ J2(^/^)Ij^5](„„ v/2)0{nv + v/2 - A), (4) 

where Jn are the integer order Bessel functions and 9 
is the step function. In this decomposition the rates 7^ 
correspond to processes in which (2n-|-l)e are transferred, 
which allows to write the dc current as Iq = e^„(2n -|- 
1)7„. In the limit v ^ 0, nv + v/2 in Eq. 10} can be 



taken as a continuous variable x. From the asymptotic 
behavior of J„ for large n one can obtain a simplified 
expression for the current at low bias 
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where cosh a = x/A. 

The predictions of Eq. © are compared with the full 
numerical results in Fig. 2. It should be stressed that the 
asymptotic expression is only valid in the voltage range 
ev < (A — A). Within this range the current is sup- 
pressed according to Eq. lO, decreasing exponentially 
with the parameter (A — A)/ev. The agreement with 
the full numerical results support our interpretation in 
terms of the AS's dynamics. 

In the absence of electron-hole symmetry (eQ 7^ 0) the 
dynamics of the AS's become more complex due to in- 
terference effects. To establish a dc current now requires 
LZ transitions between the lower and the upper state in 
addition to the mechanism previously discussed. In con- 
trast to the case of a quantum point contact (in which the 
AS's merge with the continuous spectrum at = 2n7r), 
in the S-QD-S case there is a finite probability of pre- 
serving the system coherence after a complete cycle. As 
a consequence interference effects between subsequent LZ 
transitions may appear. The two main paths for this in- 
terference correspond to the evolution following the lower 
and upper AS, as schematically shown in Fig. 1 b. The 
phase difference between these two paths is thus simply 
given by J^^ d(f>ujs{4>)/ev. One would therefore expect 

oscillations of the form cos^(2A/ew) to be observable in 
the current in addition to the ones associated with the 
onset of new MAR processes. 
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FIG. 3: Evolution of the oscillations pattern in the IV char- 
acteristics as a function of for r*/A = 5 and tg/A = 0, 
0.5, 1.0, 1.5 and 2 from top to bottom. Inset: r*/A = 2 and 
4/ A = 0, 0.2, 0.4, 0.6 and 0.8. 

Fig. 3 shows the numerical results for the low bias 
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dc current in a non-symmetric situation. As can be ob- 
served, when increasing Cg there is an overall suppression 
of the current together with the onset of an oscillatory be- 
havior of increasing amplitude. These curves correspond 
to a large r*/A value in order to disentangle these os- 
cillations from the MAR induced ones. The overall sup- 
pression is governed by the LZ probability p 0. The 
superimposed oscillations exhibit the expected behavior 
associated with the interference between AS's due to mul- 
tiple LZ transitions. A simple argument suggests that the 
amplitude of these oscillations should scale asp{l—p) and 
therefore be maximum for p =1/2 (i.e. eu/A ~ 7rr/ln2), 
as confirmed by the numerical results. When r*/A is re- 
duced one can observe a superposition of the oscillations 
arising from the two mechanisms discussed above. The 
inset of Fig. 3 illustrates the transition from one type of 
oscillations to the other in a reduced voltage range. 




FIG. 4: Zero-bias conductance defined by introducing a finite 
voltage resolution eSv/A = 0.025. The normal transmission 
from bottom to top is 0.978, 0.980, 0.982, 0.984, 0.986, 0.99 
and 1.0. 

The results presented above indicate that the IV 
curves of a S-QD-S in the Kondo regime are extremely 
non-linear in the low bias limit and thus the system can- 
not be characterized just by its zero-bias conductance. In 
an actual experiment the observability of the predicted 
behavior will be limited by several possible factors like 
the multilevel structure of the QD, environmental effects 
and inelastic relaxation processes which could be of im- 
portance at finite temperature. In the present calcula- 
tions we have assumed sufficiently small temperatures so 
that the inelastic relaxation is negligible within the volt- 
age range considered. In order to attempt a comparison 
with the results of Ref. one should take into account 
that the actual voltage resolution is limited to some mini- 
mum value Sv in such a way that the measured zero-bias 
conductance is Gs — Io{Sv)/6v. In these experiments 
6v can be estimated to be around 2.5/LtF, which corre- 
sponds to eSv/A ~ 0.025 The zero-bias conduc- 
tance defined by this finite voltage resolution is plotted 
in Fig. 4 as a function of T^/A. When Tk/A » 1 
the conductance Gs tends to saturate to a value which 



can be several times larger than the normal conductance. 
The actual value depends on 6v and is very sensitive to 
slight deviations from perfect transmission (deviations as 
small as 2% from r = 1 produce variations by a factor 
~ 10 in Gs/Gn)- On the other hand, the conductance 
exhibits an exponential suppression for Tk/ ^ < 1- AH 
these predictions are in qualitative agreement with the 
results presented in Fig. 4 of Ref. 

In conclusion we have shown that the transport proper- 
ties of a S-QD-S in the Kondo regime can be understood 
from the dynamics of the Andreev states under an ap- 
plied bias voltage. The quasiparticle current at low bias 
is controlled by inelastic transitions between the AS's 
and the extended states in the continuous spectrum. On 
the other hand, we predict novel interference effects due 
to multiple Landau-Zcncr transitions between Andreev 
states which could be worth to explore in future experi- 
ments. 
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